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Abstract 

We investigate a system of two coupled harmonic oscillators on the non-commutative 
plane M 2 . by requiring that the spatial coordinates do not commute. We show that the sys- 
tem can be diagonalized by a suitable transformation, i.e. a rotation with a mixing angle 
a. The obtained eigenstates as well as the eigenvalues depend on the non-commutativity 
parameter 9. Focusing on the ground state wave function before the transformation, we 
calculate the density matrix po(9) and find that its traces Tr(po(#)) and Tr(pg(#)) are 
not affected by the non-commutativity. Evaluating the Wigner function on M 2 . confirms 
this. The uncertainty relation is explicitly determined and found to depend on 9. For 
small values of 9, the relation is shifted by a 9 2 term, which can be interpreted as a quan- 
tum correction. The calculated entropy does not change with respect to the normal case. 
We consider the limits a = and a = ^. In first case, by identifying 9 to the squared 
magnetic length, one can recover basic features of the Hall system. 
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1 Introduction 



Because of its mathematical simplicity, the harmonic oscillator provides solvable models in 
many branches of physics. It often gives a clear illustration of an abstract idea. For instance 
a charged particle on the plane x = (xi,^) in presence of a strong uniform magnetic field 
described by the Lagrangian 

777 • 2 -* 

L = —r + A ■ x — V(x) (1) 
where the vector potential A = {A\, A 2 ) and the confining potential are 

Ai = ^BtijXj, V(x) = ^ux 2 (2) 

with 6i2 = — £2i = 1 and otherwise, gives a nice natural non-commutative system pQ and also 
a good starting point to discuss the quantum Hall effect [2|. Now it is natural to ask what 
happens if the external potential is not parabolic? This question was answered by Kim et al. 
[3 — 11] starting more than twenty years ago. They considered two coupled harmonic oscillators 
with the potential 

V(x) = - (cix\ + c 2 x\ + c 3 xix 2 ) (3) 

with constant Ci, 02,03 and explicitly determined the corresponding density matrix and the 
Wigner function as well as other quantities. 

There are many physical models based on coupled harmonic oscillators, such as the Lee 
model in quantum field theory [T2], the Bogoliubov transformation in superconductivity [13J, 
two- mode squeezed states of light [ElEJE3i the covariant harmonic oscillator model for the 
parton picture j^j, and models in molecular physics jT^]. There are also models of current 
interest in which one of the variables is not observed, including thermo-field dynamics ^7j, two- 
mode squeezed states |18| H§]. the hadronic temperature [S], and the Barnet-Phoenix version of 
information theory j20] . In all of these cases, the mixing angle a of the employed transformation 
is f, and in this situation the mathematics becomes simple. 

In this paper we study quantum mechanically a system of two coupled harmonic oscillators 
on the non-commutative (NC) plane IRj. This can be done by demanding that the spatial 
coordinates do not commute. We use the star-product to write the NC Hamiltonian and solve 
the eigenequations to get the eigenstates as well as the energy spectrum. Focusing on the ground 
state, we evaluate the corresponding density matrix. It depends on the non-commutativity 
parameter 9 but its traces as well as the entropy are not affected by the non-commutativity. Also 
we calculate the Wigner function to confirm the ^-independence of the traces. The uncertainty 
relation on Rg is found to depend on 9 and it coincides with the normal case in the limit 9 = 0. 
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For small 6, we show that this relation contains a quantum correction, which is a shift with a 
8 2 term. Also we discuss some limits of the mixing angle, namely a = | and a = 0. In the last 
case, the system can be linked to the Hall electron. 

In section 2, we give the energy spectrum and the eigenstates of a Hamiltonian describing 
two coupled harmonic oscillators. This serves as a guide in section 3 where we consider the same 
system but require that the spatial coordinates do not commute. In section 4, we deal with 
the corresponding density matrix and evaluate its traces. In section 5 we calculate the Wigner 
function on M 2 , and investigate its link with the density matrix. We determine explicitly the 
uncertainty relation as well as the entropy in section 6. Some particular cases will be considered 
in section 7. Finally, we conclude our work in the last section. 



2 Coupled harmonic oscillators 

Let us consider a system of two coupled harmonic oscillators parameterized by the coordinates 
Xi,X 2 and masses m 1 ,m 2 . This can be described by a Hamiltonian as the sum of free and 
interacting parts [TT] 

H x = 7T— p i + 7r— p i + \ {CiXl + C 2 X 2 2 + C 3 X X X 2 ) (4) 
2mi 2m 2 2 ' 

where Ci,C 2 ,C 3 are constant parameters. After rescaling the position variables 

Xi = (m 1 /m 2 )^X ll x 2 = (m 2 /mi)?I 2 (5) 

as well as the momentum 

Pi = (m 2 /mi)^i, P2 = (m 1 /m 2 )3P 2 (6) 

H\ can be written as 

H 2 = 2m ^ + P ^ + 2 + C2X ^ + C3XlX2 ) ( 7 ) 
where the parameters are 

m = (mim 2 ) 1/2 , c x = CW — , c x = CW— , c 3 = C 3 . (8) 

y mi y m 2 

As the Hamiltonian (J2|) involves an interacting term, a straightforward investigation of the 
basic features of the system is not easy. Nevertheless, we can simplify this situation by a 
transformation to new phase space variables 

y a = M ab x b , q a = M ab p b (9) 
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where the matrix 

/ cos % — sin % \ 

(M a b) =[ . I J) (10) 

\ sin 7j cos 2 / 

is a unitary rotation with the mixing angle a. Inserting the mapping (JSJ) into (J7J), one realizes 
that a should satisfy the condition 

tana = — (11) 

c 2 - Ci 



to get a factorizing Hamiltonian 

2m 



Hs = ^- (3? + ? 2 2 ) + | (e 2 "y? + e- 2 ^ 2 ) (12) 



where 



ir / 777 « ci + c 2 + v 7 (ci - c 2 ) 2 + c| 

^ = y cic 2 - c|/4, e"= (13) 



and the condition 4cic 2 > c\ must be fulfilled. 

It is convenient to separate ((H)) into two commuting parts 



# 3 = e"Wi + e-"W 2 (14) 



where and 7i 2 are given by 



Wi = ^-e^ + f e"j& H 2 = + Z-e-Oy*. (15) 

2m 2 2m 2 

First, one can see that the decoupled Hamiltonian 

H " = L qt + 1* + + V s (16) 

is obtained for rj — 0. Second, it is interesting to note that ([TflJ) can be derived by a canonical 
transformation only from 

n = n 1 + n 2 (17) 

rather than from (|14j) . This suggests that it might be advantageous to consider (JT7j) instead 
of ()14j) . Because of that Kim et al. PHIE] were focusing on the Hamiltonian ([17)1 . 

It is clear that H is a Hamiltonian of two decoupled harmonic oscillators. Thus it can simply 
be diagonalized by defining a set of creation and annihilation operators 



K t K v i _ 2 

e22/j + ^=e 2gj, o| = \/ -^—^Vi = e 2 ^ 



2hw V2mhu; 1 V 2fta; ^2mhu 

with frequency 



w = v-- ( 19 ) 

V m 



They satisfy the commutation relations 

[a i ,a]] = 5 ij . (20) 
Other commutators vanish. Now we can map Tt in terms of <2j and a\ as 

H = hw \ a\at + a\a 2 + l) • (21) 
To obtain the eigenstates and the eigenvalues, one solves the eigenequations 

H\n h n 2 ) = £m,n 2 n 2) (22) 

getting 

Kn 2 ) = Mp^|0,0) (23) 
where the ground state wave function is 

MV) = <l/i,l&|0,0) = Y^expj-^ + (24) 
as well as the energy spectrum 

£n u n 2 = { n l + n 2 + 1) • (25) 

The above solutions can be used to deduce those corresponding to H 3 , in particular its 
energy spectrum 

E 3>num = !iw (e* (ni + - V e"" (n 2 + D) (26) 
and the ground state wave function 

V>o0?) = (xi,x 2 \0,0) 



2 V 2 



^ exp 1 ~lih 



e 1 [ xi cos — — x 2 sin — J + e ' sin — + x 2 cos — 



.(27) 



While (j2I| is separable in terms of the variables y\ and y 2 , this is not the case for (|27j) in terms 
of xi and 0:2. 

In what follows, we generalize the present system to the NC case by deforming the spatial 
configuration. This will be used to investigate some physical quantities of the system, i.e. the 
density matrix and the Wigner function as well as other quantities corresponding to the ground 
state wave function (J27|) in the NC case. 
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3 Non- commutative system 



We proceed by using the NC geometry [2.1 3 to study two coupled harmonic oscillators. We 
demand that the coordinates of the plane do not commute 

[xi,Xj] = iOij (28) 

where 9ij = e^9 is the non-commutativity parameter. This relation can be obtained using the 
star-product definition 

f(x) -k g(x) = exp \ -/>,/),. <)„ \ f{x)g(y) (29) 



x=y 



where / and g are two arbitrary functions, supposed to be infinitely different iable. In what 
follows, we will use the standard commutation relations 

[Pi,Xj] = -iSij, \puPj] = (30) 



supplemented by the relation (|28j) . Together they define a generalized quantum mechanics, 
which leads to the standard one for 9 = 0. 

Now let us define the Hamiltonian (J7j) onR^. Noting that H 2 acts on an arbitrary function 
*(f,t) as 

H 2 **{r,t) = H^(f,t) (31) 
and applying the definition P9*)l we obtain 

2 / n \ 2 



9 

X2 + 2h Pl 



(32) 

With (|32p. we actually have two possibilities to get the NC version of ()17jl. This can be done 
either by transforming H 2 C via (jOJ) to obtain 

1 K K0 

n ac = — (e~\! + e\ 2 2 ) + - (e"y? + e^yl) + — (e^y 2qi - e* Vl q 2 ) (33) 



or by starting straightforwardly from ([T7jl . using (}29|) to end up with (J33|) . The effective mass 
is given by 

m 

It is useful to write (J3l?j) in the compact form 



M = (34) 



n' nc = j M {Ql + Ql) + j (n 2 + Yi) + ^ cnox - ng 2 ) (35) 
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by rescaling the variables to new coordinates Y\ and Qi. Comparing Ti.' nc to Tt, we note that 
7i' nc contains an additional term proportional to 9 which is basically the angular momentum. 
Also, the NC system shows an effective mass M which coincides with the mass m for 9 = 0. 

For the diagonalization of 7i' nc we express the position and momentum variables in terms 
of creation and annihilation operators 



which commute and satisfy the relations 



where the effective frequency 



M 

depends on 9. With the help of another set of operators 

Bl = -j=(b 1 + ib 2 ), B \ = ±={b\-ib\) 



which satisfy 



one can write 



with frequencies 



and obtain 

1 91.1 T).n f)\ = — 

\fni\n 2 \ 



(36) 



M]] = <% (37) 



, K 

Q =\T7 



B 2 = j= 2 (-b 1 + zb 2 ), Bt = ±(-b\-ibl) (39) 



[B i ,B]] = 5 ij , (40) 

n' nc = MiiB\Bi + m 2 B\B 2 + m (41) 



K9 K9 
tt 1 = tt + — ^, tt 2 = tt - -4-. (42) 



With (jUJ), we can easily solve the eigenequations 

H'* c \n l ,n 2l 9)=S'™ n2 \n 1 ,n 2l 9} (43) 



n 1 ,n 2 ,9)= i -^0f^\O,O,9) (44) 



and the eigenvalues 

S'™ na = htt ini + htt 2 n 2 + htt. (45) 
Projecting |0, 0, 9) on the plane (¥]_, Y 2 ) we find the ground state wave function 



M y,o ) = J^ P {-^(y; + y ?) }. (46, 



In terms of the x representation, it can be written as 



MQ f .I/O 
= — exp -_ 



... a .ay _„/ .a a 

e ' ( x\ cos — — X2 sin — I + e ' I x\ srn — + x<i cos — 



(47) 

We note that the analysis in the previous section is recovered for 8 = 0. 

4 Density matrix of the NC system 

Because of its relevance in determining several thermodynamic quantities, it is interesting to 
calculate the density matrix 

p(x,x') =ip(x)xp*(o?) (48) 

of the NC system. 

For the ground state wave function (|4*7)l the density matrix 

Po (x,x',d) = Mx,0)ro(?,d) (49) 

can be represented as an integral 

Po(x,x,6)= J po(x, x", 9)po(x", x 1 , 9) dx'[dx" 2 . (50) 
Tracing po(x, x', 6) over the variable x 2 , the resulting density is 

p (x 1 ,x' 1 ,6) = / tp (x 1 ,X2,0)ifj* Q (x , 1 ,x 2 ,0) dx 2 . (51) 



Evaluating this integral, we get 



poixx.x^O) = \ —— ) exp < — — (xi + Xi) sin asinh 77 



oxi> { ~ [x\ + x[ 2 ) (J cos 2 | + e-> sin 2 | ) j. ( 52 ) 



with the abbreviation 



7 = sin 2 - + e- T > cos 2 -. (53) 
'22 v ; 



The diagonal elements are 

'OM\ 1/2 f OM 



We can use f!52j) to show that the relation 

Tr(poW) = l (55) 



is satisfied as it should be for a normalized state, where the notion po(9) = po(x,x',9) is used. 
The trace of Pq{9) can be obtained by evaluating the integral 

Tr (pl(6)) = J po(x!, x[, 9)p (x[, x u 9) dx' 1 dx 1 (56) 

to end up with 

Tr(pg(0)) = (l + sinh 2 r7sin 2 a)^ . (57) 

Clearly, the traces are ^-independent and therefore they are not affected by the non-commutativity 
(|28p. The density matrix obtained by Kim et al. ^T] can be recovered by taking the limit 6 = 0. 



5 Wigner function for the NC system 

Next, we use the wave function (]37)l to determine the corresponding Wigner function, which in 
general is defined by [23] 

d 



W n (x;p) 



1 

irH 



e~^ip* n {x - s) ip n {x + s) d d s 



(58) 



for the eigenfunctions ip n (x). This definition is based on the operator formulation of quantum 
mechanics. There is another definition that can be used by introducing the ft-star-product. 
These two definitions are equivalent [21]. One concrete example was given by Dayi and Kel- 
leyane [23] who evaluated the Wigner function for an electron on the NC plane in the presence 
of a magnetic field. 

For ([37]) we evaluate the integral 

d n 

W„{x; p,6) = ( I / ( f - ? . 0) *> (2 + <>) 



to get 

W Q (x-p,9) = -f (x,9)g (p,H) 

IT 

where the functions fo(x, 9) and go(p, 9) are 

MVL 



1 



(59) 
(60) 



fo(x,9) 



9o(p,9) 



Mtt 

7th 



cxp 



exp 



h 
MVL 



,, i ct . a\ 2 __ / . a a^ 2 

e ' \x\ cos xo sin — + e ' \X\ sin — \- xo cos — 

' 2 2/ V 2 2 



h 



a 



Pi cos — — p 2 sin —j + e ' [ pi sin — + p 2 cos — 



a 



a 



(61) 



Integrating (|5U]) over the variables x 2 and p 2 , we obtain 



W (x 1 ,p 1 , 



exp 



W (x;p, 9) dx 2 dp 2 



Hy/l + sinh 2 rj sin 2 



7T 



a: 



ft(cosh r\ — sinh 77 cos a) 1 
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x 1 > exp 



ft(cosh 77 + sinh 77 cos a) 



(62) 



Of course at 9 = we recover the standard Wigner function for this system [TTj . 
One can use Wq(xi,Pi,9) to verify (J55|) and (|57ji by evaluating [9^ 



Tr(p (#)) = / W o (e)(ii,pi,0) dxxdp! 

as well as 



Tr (p 2 (6)) =2nj W$(x 1}Pu 9) d Xl d Pl . 
Indeed, these integrals confirm that the traces are independent. 

6 Uncertainty relation and entropy on Rq 

The expectation value of any operator A is defined by the matrix element 

(A) = 

where is a normalized state. Also it can be expressed as the trace 

(A) = Tr(Ap) 

over the physical states. 

The uncertainty of the operator A is given by 



AA = ^(A 2 ) - (A) 2 . 

Using either (jHSj) or (fHHjl . we calculate the product Ax\Api for the present system on 



A Xl A Pl _k / (Mtt9\_h / a 2 
'l + sinh 2 rj sin 2 a~^T + J ~~ + TTa~ 2 



where we have set a = This relation can also be obtained by using the Wigner 



mu}0 

s iicLve set u — 

tion (pD. 

For 9 = 0, we recover the standard result for the present system ^T] 



a/1 + sinh 2 rj sin 2 a 2 

For 6* — > oo, we obtain 



V 1 + sinh rj sin 2 a 
For small values of 9, (168)1 can be expanded as 

AxiAp 1 h ha 2 



\/l + sinh 2 rj sin 2 a 2 4 



F(x,p) 




Figure 1: Variation of the uncertainty product F(x, p) = 2AxiApi/h in terms of a = for 
r] = or a = 0. 

This shift can be interpreted as a quantum correction to (|69|) . In Figure 1 we show how this 
relation depends on 9. 

Finally, we calculate the entropy [22] 

S = -Tr(p\np) . (72) 

for the NC system. Using the matrix elements (po(9)) mm of po(x\, x[, 6) we determine 

S = - ^2(po(9)) mm \ia(p (9)) mm = 2 [cosh 2 // In (cosh rj) - sinh 2 ^ln(sinh?7)] (73) 

in 

which is nothing but that obtained by Kim et al. ill.. Thus we conclude that S is not changed 
by the non-commutativity ([28)1 . 

7 Particular values of a 

To discuss some limits we distinguish two different values of the mixing angle a, i.e. and 
|. In the former case we identify the non-commutativity parameter 9 to the squared magnetic 
length l 2 B . In this case we show that our system can be linked to an electron on the plane in 
the presence of the magnetic field B. 

7.1 Case a = | 

As we mentioned in our introduction, it is relevant to consider a — |, because this is the case 
for most physical systems described by two coupled harmonic oscillators. In this limit, we get 
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the ground state wave function 



the density matrix 
p (xi,x[,t 



exp ^ — \[x\ + x 2 ) cosh?] — 2xi£ 2 sinh 77 1 



7T 



h 



1/2 



nhcoshr/ 
and the Wigner function 

w {x,p,e) 



exp 



2h 



(74) 



X\ + x[)' 2 sinh^tanhr] + (x\ + xf 2 ) cosh 77 



1 



/o(a?,0) (/ O (P,0) 



(75) 



(76) 



where fo(x, 



and g (p,9) 



are 



/o(*,0) 



Mil f .\/<> r , 2 2N , , i 

exp <{ — + x 2 ) cosily — 2x 1 x 2 sinh 77 J 



71"/?, 



h 

f Mfi r . 2 ,. -1 
exp < — \[p l + p 2 ) cosh 77 + 2pip2 sinh 77 J 



The uncertainty relation simplifies to 



(77) 



y/l + sinh 2 



R /, a 2 

2V 1 + TT^' 



(7? 



7.2 Case a = 



This case corresponds to the limit C3 = 0. In the standard geometry, the system becomes 
decoupled. However in the NC plane we still have an effective coupling, which is ^-dependent. 
Before going on, we note that solving (J 13)) we have two possibilities 



,</i 



Ci 



o'l2 



c 2 



(79) 



The expressions (|79)l can be simplified by linking the present system to the Hall electron, setting 

Be 



to 



m 2mc' 



0) 



where Ib 



■g| is the magnetic length. Then 



ci = c 2 = ujy/m, e Vl = e m = 1. 



(81) 



11 



In this case, we recover the ground state wave function of the Hall electron 

ipo(x, Oh) 



a=0 



mu 



This leads to the density matrix 

Poixx^'^On) 



cxp 



(^) 1/2 exp{^ 1 +*?)}. 



mu 



2h 



2 , 2^ 
Xi -+- X 2 , 



a=0 



The corresponding Wigner function reads 



w (x,p,e H ) 



= ~fo(x,0 H ) go(p,0 H ) 

a=0 7T a=0 



a=0 



(82) 



53) 



where f (x, H ) and g (p, H ) 

a=0 



arc 



o=0 



fo(x, H 
9o(p,0 H ) 



a=0 



a=0 



mu r mu 

— exp rir 

mu 



mu 



cxp 



h 



(Pl + Pl)}- 



55) 



8 Conclusion 



We have investigated quantum mechanically a system of two coupled harmonic oscillators on 
the non-commutative plane by requiring that the spatial coordinates do not commute and 
employing the star-product definition. By writing down its NC Hamiltonian and making use 
of a suitable transformation, i.e. a unitary rotation with the mixing angle a, we have ended 
up with a diagonalized system where the condition (jllj) for a was taken into account. By 
solving the eigenequations, the eigenstates and the energy spectrum as well as the ground state 
wave function are found to depend on the non-commutativity parameter and to coincide for 
= with those for the standard case. This was used to determine explicitly the ground state 
wave function of the NC system before the transformation, which was the starting point of our 
interest. 

Subsequently, we have used the above tools to determine some physical quantities corre- 
sponding to the ground state. We have evaluated the corresponding density matrix, which is 
^-dependent. However, we have shown that its traces, i.e. Tr(p (#)) and Tr (pq(6>)), are not 
affected by the non-commutativity. Also the Wigner function was calculated for the NC system 
and used to confirm the trace properties. We have explicitly evaluated the uncertainty product 
and some relevant limits were investigated. By taking small values of 0, we have found that 
the standard relation was shifted by a term proportional to O 2 . This effect was interpreted as 
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quantum correction to the normal case. Of course this shift vanishes for 6 = 0. We have found 
no dependence for the entropy of the NC system, it is the same as that for the normal case. 

Finally, we have considered some limits of the parameters. In particular, we have studied 
the case a = | which is relevant for many physical systems. In another interesting case, a = 0, 
the standard system becomes decoupled. However, the NC system retains an effective coupling 
due to the parameter 9. By identifying 6 to the squared magnetic length, basic features of the 
Hall system were recovered, in particular its ground state and its Wigner function were given. 

This work should be a good starting point to investigate different issues related to the 
physical systems described by two coupled harmonic oscillators. In particular, it can be used 
to investigate different models proposed to study many physical problems as mentioned in our 
introduction, see [30111 and [12 — 20], by deforming the spatial coordinates. We will return to 
these issues and related matter in future. 
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